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SYMBOLS 

C sum of the observed stimulus levels in first n blocks 

n (see Eqs. 2 and 3) 

d step size (distance between adjacent stimulus levels) 

E efficiency (see Eq. 9). 

F(x) probability of correct response due to sensory mechanism 
(see Eq. 4) 

G(x) probability of correct response due to guessing (see Eq. U) 
k trial block length 

n number of trial blocks 

P(x) probability of correct response (see Eq. 5) 

P_ probability of decreasing the stimulus level following a 

block of trials 

P + probability of increasing the stimulus level following a 

block of trials 

P q probability of remaining at the same stimulus level following 

a block of trials 

r decrease stimulus level if number of correct responses in a 

block is > r 

RE relative efficiency isee Eq. 10) 

s increase stimulus level if number of correct responses in a 

block is < s 

s a general fixed stimulus level 

x^ ith stimulus level fixed prior to an experiment 

y^ stimulus level used in (j+l)st trial block 

y initial stimulus level 

o 

<j> phasing factor; distance from y to the fixed stimulus level 

closest to y 

y true threshold value 
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0 estimator of y (Eg. l) 

asymptotic frequency of stimulus level i 
a standard deviation of the underlying normal distribution F(x) 

SUMMARY 

The block up-and-down, two-alternative, forced-choice experi- 
mental design for the estimation of sensory thresholds is investigated. 
A mathematical model of the procedure is developed and recursion for- 
mulas are derived for bias and mean-square error of the sample average 
estimator. Block designs for various step sizes are compared on the 
basis of two measures of efficiency: (l) efficiency expressed as the 

reciprocal of the mean-square error of the average estimator per trial 
and (2) relative efficiency of the average estimator with respect to 
the probit estimator based on fixed stimulus levels. 

INTRODUCTION 

Up-and-down or "staircase" designs have been extensively analyzed 
and applied during the past few years to psychophysical experiments. 1-8 
These experimental designs, widely used in other fields such as bio- 
assay and explosives research, 9 " 13 were developed to estimate points 
on a response function where responses are quantal; that is, responses 
are categorized as occurring or not occurring following a stimulus. 

In the application considered in this report, this function is the 
probability of correct response and is assumed to be in the form of a 
gaussian psychometric function adjusted by a correction for guessing. 
The up-and-down method consists in giving the experimental subject a 



series of stimuli at a fixed number of sequential trials according to 
the following rule: (a) following a nonresponse, increase the stimu- 

lus to the next higher stimulus level for the next trial, (b) follow- 
ing a response, decrease the stimulus to the next lower stimulus level. 
This procedure tends to concentrate testing around a fixed point on 
the response function whose corresponding stimulus level may be defined 
as a "threshold," or point of subjective equality. The block up-and- 
down method is an extension of the classical up-and-down method in 
which the decision to raise or lower the stimulus level is based on 
the outcome of a block of several trials rather than on the outcome of 
just one trial. 

Recent work in this area has been concentrated on developing 
efficient sequential methods - experimental designs in which the num- 
ber of trials is a random variable that depends on the past history of 
the trial sequence . 6 ’ 9 5 12 5 14 9 1 5 Nonsequential designs are sometimes 
easier to apply experimentally, however, and are still used in 
psychometric testing . 1 ’ 4 ’ 5 

Up-and-down designs , when used to determine sensory thresholds 
for human subjects, are usually applied within the framework of a 
two-alternative, forced-choice procedure; that is, the Subject is 
required to give one of two types of response to a stimulus presented 
in the two categories with equal probabilities. For example, Clark 
and Stewart , 5 in a study of angular acceleration, required a subject 
in a centrifuge to respond right or left, depending on which direc- 
tion he subjectively perceived that he was accelerating. The 
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difficulty in using an up-and-down design with the restriction of 
forced choice is that the probability of responding "correctly" to 
a stimulus varies from 0.5 to 1 as the stimulus level is increased, 
rather than from 0 to 1 which is assumed in other contexts, such as 
bioassay. 10 ’ 12 ’ 13 The probability of 0.5 at zero stimulus level 
arises from the assumption of random guessing under the restriction 
of an equally probable forced (binary) choice. 

Experience with empirical response curves has shown that this 
restriction results in a probability-of-correct-re^ponse function 
that is not symmetric about the threshold; that is, the probability 
of increasing a step when below threshold is not equal to the proba- 
bility of decreasing a step when above threshold. This lack of sym- 
metry can also be deduced from the mathematical model to be developed 
below. Since classical up-and-down designs assume symmetry of the 
response function, we undertook to determine how the asymmetry induced 
by the forced-choice technique affects bias and precision of the 
average estimator of threshold. This method involves averaging 
observed stimulus levels over a series of trials 10 ’ 13 to estimate 
threshold, and is found quite frequently in the psychometric 
literature . 1 > 4 > 5 > 7 > 8 

This report summarizes an investigation of the bias and precision 
of the average stimulus level estimator (the dose average estimator 
of Ref. 13) for a block up-and-down, two-alternative, forced-choice 
design (BUDTIF). A mathematical model of the forced-choice procedure 
is developed, and well-known recursion formulas for exact bias and 
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mean-square error of the average estimator are modified to handle 
the forced-choice case. Since this was not intended to be an exten- 
sive parametric study, the main results are presented in terms of 
optimal block design, that is, the number of trials at a given stimu- 
lus level and the appropriate decision proqedure to raise or lower 
the following stimulus level, which maximizes some measure of effi- 
ciency of the average estimator. The measure of efficiency chosen, 
and one that incorporates the influence of most of the relevant 
parameters, is the reciprocal of the me an- square error of the average 
estimator per trial. 

The parameters that influence the up-and-down procedure are: 

(1) initial stimulus level, (2) step size (i.e., the fixed distance 
between stimulus levels), (3) number of trial blocks, ( k ) block 
design (described below) , and (5) phasing factor. All these parameters, 
except the phasing factor, will be discussed. The phasing factor, 
which is the distance from the threshold to the stimulus level nearest 
the threshold, was found to have a slight but insignificant effect on 
mean-square error and bias for moderate trial sequence lengths. 

Because of this small effect on mean-square error and bias of the 
threshold estimator, and since the phasing factor cannot be known to 
the experimenter in practice, all results are presented for a phasing 


factor of zero. 
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I. MODEL OF THE EXPERIMENTAL PROCEDURE 

The general psychophysical method of interest is an expansion 
of the BUDTIF procedure developed by Campbell 8 and the multiple up- 
and-down (MUD) procedure used in bioassay. 1 1 5 12 The rules for 
manipulating the independent variables are: 

1. Choose a set of stimulus levels that are equally spaced 
(usually in log units of physical magnitudes). 

2. Perform a sequence of trials in blocks of length k. After 
each block of trials is completed at a given stimulus level, 
select the stimulus level for the next block as follows: 

a. Increase to the next higher level following s or fewer 
correct responses in the present block. 

b. Decrease to the next lower level following r or more 
correct responses (r > s). 

c. Remain at the same stimulus level following a number of 
correct responses between s and r (not including s or r) . 

3. Terminate the experiment after n blocks of trials. 

These parameters determine the block design (k, s, r). 

This rule for changing the stimulus level after a trial block is 
more general than Campbell's in that he considered only the case in 
which the decision to raise, lower, or keep the stimulus level the same 
was based on whether the proportion of correct responses (x/k) was less 
than, greater than, or equal to the desired proportion (p) of correct 
responses, respectively. In confining himself to this rule, he assumed 
that the condition x/k = p had to be a possible outcome, where x is 



8 


the number of correct responses in a block of k trials , and p is 
some target percentage correct which is being tracked. This .con- 
dition necessarily confined him to block lengths that are multiples 
of 1+ when p = 0.75, an unnecessary restriction when the above rule 
is used. 

Various suggestions have been made concerning the best way to 
avoid the bias that results in the threshold estimate when the first 
stimulus value is far away from the true threshold. Campbell 8 sug- 
gested using only those levels, in the threshold calculation , that 
have been used at least twice within the trial sequence ; Brownlee 1 3 
suggested calculating threshold from those levels used after the first 
reversal (change of direction) of the staircase; Hsi 12 suggested 
beginning a sequence with single-trial blocks and then switching to 
k-trial blocks after the first stimulus reversal. Since all these 
methods are designed to place the initial trial block for the cal- 
culation of the threshold estimate in the vicinity of the threshold, 
a suggestion of Hsi 12 was followed and in' the computer study the 
starting stimulus level was confined to within three standard devia- 
tions of the true (simulated) threshold. 

II . STATISTICAL METHOD 

Recursion Formulas for Bias and Mean Square Error 
The average estimator is defined as follows: Let be the 

ith stimulus level, fixed prior to the experiment, i = 1, 2, . . ., 
where x^ +l > x i for all i, and let y , j = 0, 1, 2, . . . , be 
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the stimulus level used in the j+1 trial block. If p is the true 
threshold, consider estimating p by p, where 

n n 

P = XI 3 2 (i) 

j=l J j=l J 

and k is the number of trials per block in a sequence of n blocks . 

As in Reference 13, the initial stimulus level y is not included 

o 

since it was chosen by the experimenter, and the level y n is 
included since it was selected on the basis of the experiment. 

Recursion formulas for the exact bias and mean-square error of 
p are given by.Hsi. 12 Hsi generalized to the case of k > 1 the 
formulas developed by Brownlee et al. 13 for k = 1. The only modifi- 
cation of Hsi's formulas that is required is the definition of P + , 

P Q , and P , the probabilities of increasing the stimulus level after 
a block of trials, staying at the same level, and decreasing the 
level, respectively. These probabilities will be derived later. The 
recursive formulas for the bias and mean-square error, as given in 
Reference 12, are as follows (given that p = 0): 

E[ c n+ i (y ) ly Q =x i ^ = E[Ci(x.)] + P + E[C n (x. + 1 )] + P^tC^x.)] (2) 

E[ c n+ i(y)|y 0 =xJ] = E 2 [C 1 (x i )]+P + E[C^(x i+1 )]+P o E[C^(x i )] 

+P E [ C 2 ( x )]+2kx,{E[C (x )]-E[Ci(x. )]} 

- ni-l i n+li l 

+2kt(x i + r x i>V E t c n ' x i + i )] - ( Wl )P - E[0 n U i-l )]1 


( 3 ) 
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where 

C n (y) = 2 

j=l J 

A computer program was written to compute the exact bias and mean- 
square error using the above formulas . On a "per trial" basis , the 
result is 

Bias = E[C n+1 (y) |y Q =x i ]/kn 
MSE = E[C n+1 (y)|y o =x.] 2 /k 2 n 2 

Statistical Model 

The probabilities P + , P Q , P may be derived as follows: In a 

forced-choice procedure, a correct response at a given stimulus level 
x is the result of a correct response due to the actual sensory mech- 
anism under study, or is due to some random response process. If the 
probabilities associated with these two events are represented by 
F(x) and G(x) , respectively,, and if we represent the probability of a 
correct response by P(x), then 

P(x) = F(x) + [1 - F(x) ]G( x) (k) 

where 0.5 £ G(x) < 1 and 0 < F(x) < 1. The psychological interpreta- 
tion of this equation is:, with probability F(x) , the subject detects 
the signal (perceives the stimulus); with probability 1 - F(x), he 
does not, and if he does not then he responds with a probability of 
G(x). The inequalities imply that 0.5 < P(x) < 1. The case con- 
sidered in this report occurs when G(x) is independent of x and is 
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completely random (i.e., guessing). In such case, G(x) =0.5 for 
all x and P(x) = [l + F(x)]/2. The term F(x) is referred to in 
the literature as the psychometric function and is commonly assumed 
to he described by the normal cumulative distribution. The same 
assumption is made here ; thus , 


P(x) 



(2ir a 2 )"" 1 / 2 exp 


r ifc-A 2 " 

L 2 \ / - 


dt \ 


(5) 


Since F(x) represents the actual sensory mechanism under study, by 
definition its mean (in this case also its median) y is the threshold 
value of interest. If y is substituted for x in Eq.. 5» 

P(y) = 0.75. This is the reason for the 75$ target percentage correct 
that is commonly estimated in forced-choice situations. 2 

Now, the probability of increasing, decreasing, or keeping the 
same stimulus level following a trial block may be represented by the 
binomial distribution partial sums. Thus, 


P + = probability of increasing the stimulus level following a block of 
length k 

= probability of s or fewer correct responses in the block 

-VQtpwm - P<*»- 

m=o 


P = probability of decreasing the stimulus level following a block of 
length k 

= probability of r or more correct responses in the block 
k 

= 23(») [PU)]m[1 ‘ p(x,lk " m 

m=r 
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P = probability of remaining at the same stimulus level following 
a block of length k 

= probability of between s and r correct responses 

= 1 - P - P 

+ — 

Without loss of generality, the standardized form is assumed (y and 
a were taken to be 0 and 1, respectively). The step size d was 
taken to be constant for a given trial sequence (d = x^ - x^ for 
all i) . 

Asymptotic Frequency Distribution of Stimulus Levels 
The up-and-down method can be thought of and modeled as a random 
walk on the real line with fixed step size. This way of treating the 
problem is useful for deriving the asymptotic frequency distribution 
of stimulus levels, that is, the relative frequency with which each 
stimulus level is visited in an infinitely long trial sequence. 
Tsutakawa 11 derived this distribution for the bioassay case, and the 
only modification of his formulas that we require is the interpreta- 
tion of P + and P , the probabilities of increasing a step and decreas- 
ing a step. As shown above, these are the probabilities of the tails 
of a binomial density function. 

If we let p* be the value of p that makes P + = P , then it 
can be shown that for p* > 1/2: 

P + > P if p < p* 

P + < P if p > p* 

indicating that the stimulus series will not drift to plus or minus 
infinity, and that the asymptotic distribution of stimulus levels will 
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have a finite mean. Therefore, Tsutakawa's formulas can he applied. 
Table I lists the asymptotic distributions that were calculated for 
different combinations of k, s, and r, along with their means and 
variances. Since y converges in probability to the asymptotic mean, 
for y = 0 this mean is the asymptotic bias of the block up-and-down 
procedure (as n -*■ 00 ). 

III. RESULTS AID DISCUSSION 

Formulas 1 and 2 were used to compute exact bias and mean-square 
error (MSE) of the threshold estimator y. The parameters that 
determine bias and precision of y are: (l) the block design (k,s,r), 

(2) the initial stimulus level y , (3) the step size d, (1+) the num- 
ber of trial blocks n, and (5) the phasing factor. The phasing factor 
was discussed earlier and eliminated from consideration. Although the 
remaining four parameters are under the control of the experimenter, 
he seldom has enough prior information about the true values of y 
and cr to allow him to select y , d, and n optimally. This study 
concentrates, therefore, on finding block designs that are good over 
wide ranges of y , d, and n. Some sample curves of bias and MSE are 
included, in special cases, to provide some insight into the relative 
influence of these parameters. 

The block design (k,s,r) will be considered first. Not all com- 
binations (k,s,r) are feasible. Since the up-and-down method concen- 
trates testing around a so-called target percentage correct, we must 
see what percentile of the function P(x) is actually tracked by the 



procedure. If the two tail probabilities of the binomial distribution 
discussed in the previous section are equated, and the resulting equa- 
tion is solved for the binomial parameter p, the result is the prob- 
ability of correct response tracked by the up-and-down procedure. 

Table II gives this value of p as a function of (k,s,r) along with 
the normalized stimulus level that yields the probability p; that is, 
the value of z is found by numerically solving the equation 



for z, given p. The stimulus level that yields p is then 
x = p + zo. Because of the asymmetry of the response function, x is 
not exactly equal to the asymptotic value of p listed in Table I. 
However, these two numbers are expected to be close for suitable step 
sizes d; combinations (k,s,r) chosen from Table II for further investiga- 
tion were those that yielded values close to 0.75 for p. As shown in the 
derivation of P(x), the inequality of 0.5 < p < 1 must hold, which 
accounts for the missing entries (k,s,r) in the table. The futility 
of considering combinations (k,s,r) that do not appear in Table II 
can also be shown by the asymptotic behavior of the up-and-down 
series in this case. It can easily be shown that for p < 0.5» where 
p is the probability that equates the two tail probabilities dis- 
cussed above, the probability of the stimulus level x decreasing on 
the next trial following a given trial is greater than the probability 
of increasing (i.e., P > P + ). Therefore, the 
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process tends to drift to and the absolute bias and MSE increase 
without bound. 

A reasonable way to choose the best s and r for a given block 
size, k, or to choose the best k for a fixed total number of trials, 
is to base comparisons on the amount of information (reciprocal of 
MSE )• per trial . 13 This measure is called the efficiency E: 



where k is the block size, n is the number of trial blocks, and 
MSE is the observed mean-square error averaged over all initial 
stimulus levels, y Q . Figures 1 through 3 give E as a function of 
the total trial sequence length nk for step sizes d = 1/2, 1, and 
2. The values of s and r chosen for each k were those that gave 
the highest efficiencies over the range of trial lengths considered. 
Whenever no single combination of s and r was uniformly best for 
all trial lengths, one of the two best was plotted. The figures show 
very strikingly the relatively poor performance of k = 2, the value 
of k very often used in psychometric work, especially for long sequence 
lengths and large step size. On the other hand, the uniformly good 
performance of k = 5 is surprising. Note that a block size of 5 
would be ruled out if the criterion stated on page 1177 of Reference 2 
were followed (the rule that it must be possible for the proportion of 
correct responses out of the block of five to be exactly equal to 0,75). 

Another way of looking at the up-and-down method is to compare 
its performance to that of a fixed stimulus level design (e.g., the 
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classical z-score method of threshold estimation). A fixed- level 
method very similar to the z-score method, and that has served as a 
basis for comparison with the up-and-down method in bioassay, 1 0 > 12 , 1 3 
is the method of probits developed by Finney. 16 This method is based 
on a maximum likelihood estimate of threshold, and is therefore known 
to be asymptotically efficient. The fixed-level design chosen was 
based on an equal number of trials at each of five stimulus levels , 
and the formula for the variance of p given in Reference 1 6 was 
used, with a assumed to be 1 and the weights adjusted to conform to 
the response function in this report, P(x). The relative efficiency 
(RE) of the BUDTIF method is defined as 

MSE (PROBIT) / s 

MSE (BUDTIF) U0; 

Representative efficiency curves in Fig. 4 show that the BUDTIF method 
is relatively efficient for starting levels away from threshold and 
for small trial sequence lengths. Even for starting levels at threshold, 
the BUDTIF method seems to be more efficient than the probit estimator 
for large block sizes. Relative efficiencies were averaged over start- 
ing levels, y o , for various (k,s,r) and are shown in Table III for 
step size d = 1. 

Bias and MSE of the estimator p are shown in Fig. 5 for k = 8 
and d = 1.5. Mean-square error based on a constant total trial 
sequence length (nk = 20) is plotted in Fig. 6 for various combina- 
tions (k,s,r). It is obvious that starting stimulus levels far away 



17 


from threshold inflate the bias and MSE; however, it is not uniformly 
true that threshold is the best starting position. Because of the 
particular response function assumed, the bias for starting levels at 
threshold tends to be negative and becomes worse as the trial sequence 
length increases. The block size of k = 8, however (Fig. 5), is 
remarkably stable and has a uniformly small positive bias up to moderate 
sequence lengths. Table IV shows that the asymptotic mean stimulus 
level is negative for most values of (k,s,r), explaining why the bias 
goes negative with increasing trial sequence lengths. 

The effect of step size d, as shown in Fig. 6, is similar to that 
previously shown for the bioassay case. If y Q is far from threshold, 
the MSE increases with d. A small d gives a very precise threshold 
estimate if is close to p, but a large d is more stable with 

respect to MSE as y Q gets farther from y. 

IV. CONCLUDING REMARKS 

Much of the past research on sampling properties of statistical 
estimates of sensitivity thresholds has been based on lengthy and 
inaccurate Monte Carlo computations that required repeated sampling 
of the responses of a simulated human subject. Although the Monte 
Carlo approach is necessary for analyzing some of the proposed 
sequential procedures that are analytically intractable, the methods 
described in this paper are adequate for computing exact bias and 
precision for fixed-length block designs in which an average estimator 
is used. This makes possible the rapid search of many more parameters 
than is practical with Monte Carlo procedures . 



18 


This investigation of the block-up-and-down design has shown that 
the design often used in the past, based on a trial sequence length of 
two, is one of the worst from the standpoint of sampling efficiency. 
Efficiency curves were provided to assist the experimenter in choosing 
an appropriate block design. 

REFERENCES 

1 E. G. Heinemann, "The Relation of Apparent Brightness to the 
Threshold for Differences in Luminance," J. Exp. Psychol. 6l(5) , 

389-399 (1961). 

2 C. ¥. Stuckey, C. L. Hutton, and R. Campbell, "Decision Rules 
in Threshold Determination," J. Acoust. Soc. Amer. 40 ( 5 ) , 117*1—1179 
( 1966 ). 

3 T. N. Cornsweet, "The Staircase Method in Psychophysics," Amer. 

J. Psychol. 75, U85-U91 (1962). 

4 S. Feeny, P. K. Kaiser, and J. P. Thomas, "An Analysis of Data 
Gathered by the Staircase Me'thod," Amer. J. Psychol. 79., 652-65** (1966). 

5 B. Clark and J. D. Stewart, "Comparison of Three Methods to 
Determine Thresholds for Perception of Angular Acceleration," Amer. 

J. Psychol. 81, 207-216 (1968). 

6 M. M. Taylor and C. D. Creelman, "PEST: Efficient Estimates on 
Probability Functions," J. Acoust. Soc. Amer. Ul(U) , 782-787 (1967). 

7 R. A. Campbell and E. Z. Lasky, "Adaptive Threshold Procedures: 
BUDTIF , " J. Acoust. Soc. Amer. 4U(2) , 537-5^1 (1968). 

8 R. A. Campbell, "Detection of a Noise Signal of Varying Duration," 
J. Acoust. Soc. Amer. 35(ll) , 1732-1737 (1963). 



19 


9 G. B. Weatherill, "Sequential Estimation of Quanta! Response 
Curves," J. Roy. Statist. Soc., Ser. B 25, 1-48 (1961). 

* 

10 W. J. Dixon and A. M. Mood, "A Method for Obtaining and Analyzing 
Sensitivity Data," J. Amer. Statist. Assoc. 43, 109-126 (1948). 

13 -R. K. Tsutakawa, "The Random Walk Design in Bio-Assay," J. Amer. 
Statist. Assoc. 62 ( 319 ) > 842-856 (1967). 

12 B. P. Hsi, "The Multiple Sample Up-and-Down Method in Bioassay," 
J. Amer. Statist. Assoc. 64(325) , 147-162 (19 69). 

13 K. A. Brownlee, J. L. Hodges, and M. Rosenblatt, "The Up-and- 
Down Method With Small Samples," J. Amer. Statist. Assoc. 48 , 

262-277 (1953). 

ltf G. B. Wetherill and H. Levitt, "Sequential Estimation of Points 
on a Psychometric Function," British J. Math. Statist. Psychol. 18 , 

1-10 (1966). 

15 M. M. Taylor and C. D. Creelman, "PEST: A Rapid Technique for 
Finding Arbitrary Points on a Psychometric Function," Psychonom. Soc., 
Chicago (1965). 

16 D. J. Finney, Pvobit Analysis (Cambridge Univ. Press, 1952). 



TABLE I. Asymptotic frequency distributions of stimulus levels for step size d = 1. Triplets 
in parentheses are (k,s,r); phasing factor <j> = 0.0. 
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TABLE II. True probability, p, of correct response 
tracked by the design (k,s,r); stimulus level that yields 
this p is z. 
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TABLE III. Relative efficiency of the average 
estimator vith respect to the probit estimator 
averaged over all starting levels y ; d = 1; 

RE = MSE (PROBIT) /MSE(BUDTIF) . 


n 


(k,s,r) 5 10 20 25 50 


( 2 , 1 , 2 ) 

4.36 

( 4 , 2 , 4 ) 

3.76 

( 5 , 3 , 5 ) 

3.56 

( 6 , 4 , 6 ) 

2.93 

( 8 , 4 , 8 ) 

2.19 

( 10 , 6 , 9 ) 

— 


3.03 

1.72 1.47 

0.94 

3.63 

3.22 


4.10 

5.03 — - 

— 

3.26 

. 

— 


3.26 
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TABLE IV. Means and variances of the asymptotic frequency 
distribution of stimulus levels . 


(k,s ,r) 

d = 1/2 

d = 

1 

d = 1-1/2 

Mean 

Var. 

Mean 

Var. 

Mean 

Var. 

(2,1,2) 

-0.3552 

0.6142 

-0.5408 

1.6412 

-0.7557 

3.1476 

(3,2,3) 

.1686 

.4210 

.1045 

1.0267 

.0320 

1.838 

(4,2,4) 

-.1306 

. 2737 

-.1938 

.6764 

-.2653 

1.254 

(5,3,5) 

.1448 

.2340 

. Ill8 : 

.5725 

.0853 

1.056 

(6,3,6) 

-.03663 

.1784 

-.07123 

.4382 

-.1132 

.8072 

(7,5,6) 

.0691 

.3081 

.02839 

.7617 

.0014 

1.399 

(8,4,8) 

.0188 

.1327 

-.005628 

.3192 

-.0368 

.5658 

(9,4,9) 

-.0862 

.1127 

-.1082 

.2550 

-.1390 

.4335 

(10,6,9) 

-.0446 

.1649 

-.07208 

.4245 

-.0943 

.8293 
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FIGURE LEGENDS 

Fig. 1. Efficiency E = l/kn (MSE) versus trial sequence length (kn) 
for step size d = 0.5. Mean-square error (MS E) was averaged 
over all starting stimulus levels y Q . Numbers in parentheses 
are (k,s,r). 

Fig. 2. Efficiency E = l/kn (MSE) versus trial sequence length (kn) 
for step size d = 1. Mean-square error (MSE) was averaged 
over all starting stimulus levels y Q . Numbers in parentheses 
are (k,s,r). 

Fig. 3. Efficiency E = l/kn (MSE) versus trial sequence length (kn) 
for step size d = 2. Mean-square error (MSE) was averaged over 
all starting stimulus levels y Q . Numbers in parentheses are 
(k,s ,r ) . 

Fig. 4. Relative efficiency (RE) of the average estimator with 

respect to the probit estimator for nk = 40. Triplets in paren- 
theses are (k,s,r) and d = 1. 

Fig. 5. (a) Bias of the estimator y versus initial stimulus y Q 

and sample size n; (k,s,r) = ( 8 » U s 8 ) and d = 1.5« (b) MSE of 
the estimator y versus initial stimulus y Q and sample size n 
(k,s,r) = (8,4,8) and d = 1.5. 

Fig. 6. Mean-square error of the average estimator y for constant 
trial sequence length (nk = 20) for k = 2, 4, and 10. Numbers 
in parentheses denote (k,s,r). 
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Figure 2 
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